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Abstract. The following nonlinear PDE: 
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derived from 2D Euler equations, is shown to admit smooth similarity solutions, which as 
t T~ create shocks of the type sign a: and H(—x) and also other blow-up singularities. 
Some of the blow-up solutions are shown to admit unique extension beyond, for t > T. 
All the similarity reductions lead to singular boundary value problems for 2nd-order 
ODEs. More complicated solutions with blow-up angular swirl are discussed. 

1. Introduction: Euler equations and a related PDE with shock waves 

1.1. Euler equations and shocks: beginning of derivation. Consider the 2D Euler 
equations (the EEs) of inviscid incompressible fluids 

(1.1) ut + {u-V)u = -Vp, divu = in x {0,T) {u={ui,U2f), 

with bounded initial divergence-free Uq of infinite kinetic L^-energy. Using Leray's for- 
mulation with the projector P = / — VA~^(V-), (11.11) is written in the form fL2}, p. 30] 



is the 2x2 Jacobian matrix of the solenoidal vector fleld u. 

It is well-known that, even in dimension 2 (to say nothing about the 3D ones), the 
standard notion of weak solutions is not sufficient for establishing the unique solvability of 
the problem in the case of the unbounded initial vorticity, where curl uq ^ L°°. See Majda- 
Bertozzi [T2l § 8] and recent surveys in Bardos-Titi [1], Constantin [4j, and Ohkitani 
[13] for further details and key references, as well as Pomeau et al [17] for discussion of 
somehow related open questions on shock waves in hyperbolic systems. It seems obvious 
that natural reasons for this difficulty are two fold. Firstly, there is no still a proper 
definition of "entropy solutions" for the EEs admitting sufficiently strong "shock waves" . 
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(1.2) 




^^tr (Vu(y,t))2dy, where 



In other words, general principles of formation of various singularities of shock wave types, 
as well as L°°-blow-up, are in general unknown. Secondly, on the other hand, even the 
corresponding system without the nonlocal term 

(1.3) Ut + {Du)u = in x M+ 

is a rather difficult divergent hyperbolic system 2x2 with two independent spatial variables 
X = {xi,X2), so it does not automatically fall into the scope of modern entropy theory 
regardless of its definite recent achievements; see Bressan [2] and Dafermos [S]. Of course, 
some classes of solutions of systems such as fll.3p can be represented via characteristics (see 
interesting physically motivated formal discussions in tl2l|), though, in general, choosing 
right weak "entropy" solutions is not straightforward and is an open problem. 

Our goal is to introduce a new PDE model, which may refiect some features of formation 
of finite-time singularities in systems such as f 1 1.1 1) (and partially (II. 3p ) and, in addition, 
having vorticity-swirl structure that can be somehow adequate to the EEs. To this end, 
to discuss possible ways of shock wave formation for the EEs, let us write down (11. ip in 
the polar coordinates (r, ip) , 

fUt + UUr + lVU^-lV^ = -Pr, 

(1.4) I Vt + UVr + lVV^ + lUV = ~'-P^, 

[Ur + lU + lV^ = 0. 

Though (11.41) is a system of three PDEs, which is equivalent to the nonlocal equation 
(II. 2p . we consider a rather hypothetical situation of formation of shocks created by local 
differential operators onljQ, though, as we have seen, this scenario is expected to be also 
relevant to the hyperbolic system (II. 3p . 

Let us then more precisely define a manifold of solutions with shocks to be studied: 

(i) Formation of bounded shocks (or bounded rarefaction waves) is not essentially af- 
fected by the pressure terms, i.e., by the nonlocal members (however, the div-free equation 
remains important); and 

(ii) As in the classic case (ll.7p below, two quadratic first-order operators in the first 
and in the second equations in (11.41) are leading together with two differential terms in 
the div-free equation (the latter is reasonable for finite shocks): 

fUt + UUr + ... = 0, 

(1.5) lVt + lVV^ + ... = 0, 

[Ur + lV^ + ... = 0. 

Thus, we omit other operators that are expected to be bounded on such singularities. 



Obviously, this is not the case for 3D unbounded L°°-singularities, where quadratic nonlocal terms 
must be involved and this settles the core of such a remarkable open problem, [1]. 
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Replacing V by rV for convenience yields 



;i.6) 



V ^rV 



'Ut + UUr + ... = 0, 
Vt + VV^ + ... = 0, 

Ur + V^ + ... = 0. 



The first two equations are indeed the famous ID Euler equations from gas dynamics, 

(1.7) ut + uu^ = {] in MxM+, u{x,Q) = uq{x) e n 1°^ , 

whose entropy theory was created by Oleinik [131112] and Kruzhkov [TT] (equations in M^) 
in the 1950-60s; see details on the history, main results, and modern developments in the 
well-known monographs [21 [5l I20j^ . Typical formation of shocks for (11. 7p as t — >^ 0~ (the 
blow-up time here is T = 0) is described by the canonical self-similar blow-up solution 

1 for X < t, 



;i.8) 



U- 



{x,t) 



1 fo^ 



l-l < 1 

X > —t. 



-1 for 

Therefore, the following steady shock is formed: 

(1.9) lim^^o- u_{x,t) = S^{x) = —signx, 

with convergence in L^(M) and uniformly on any subset bounded away from the origin. 



Moreover, we have that, as t — 



(1.10) 



DxU- 



{x,t) 



-oo for any |a;| < \t\ and 



DxU. 



.{x,t) -25{x). 



Indeed, such a standard scenario of shock wave formation is not applicable for the 
system in (11. 6p . since the last div-free equation 

(1.11) Ur + V^ + ... = Q 

prohibits simultaneous blow-up of derivatives Ur and to — oo. Therefore, the system 
(II. 6p needs another more involved treatment of possible (if any) shock waves. We next 
introduce such a nonlinear PDE model. 



1.2. A simplified PDE model and main results. 

tions in (11.61) on the "critical manifold" (11.111) . i.e.. 



Thus, we consider first two equa- 



;i.i2) 



Ut + UUr 

Vt + vv^ 



0, 










= VUr, 



with given bounded smooth initial data {Uq, Vo)(r, Lp). Writing this in the form 

{\n\U\)t = V^, 
X\n\V\)t = Ur, 



;i.i3) 



^First study of finite time singularities as shock waves of quasilinear equations was performed by 
Riemann in 1858 [18] (Riemann's method and invariants are originated therein); see [H [17] for details. 
The implicit solution of the problem ()1.7|) . u(x, t) — uq{x—u{x, t)t) (containing the key wave "overturning" 
effect) was obtained earlier by Poisson in 1808 [16j; see [T^- 
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we find on integration of tlie second equation tliat 



(1.14) in|X|)| =/J/7,(s) =^ V{t) = VoeIoUris)_ 

Substituting tliis V{t) into tlie first equation for U yields tlie following nonlocal PDE: 

(1.15) Ut = UVo^elo Ur(s) + (e/o ^'■W) /J [/,^(s) = J, + J^. 

According to (11. 150 . there exist at least two different mechanisms of shock wave formation 
corresponding to nonlinear integral operators Ji and J2, and possibly, other types via their 
nonlinear interactions. Note that the problem corresponding to J2 is assumed to be truly 
2D (and hence difficult), while that for Ji can be reduced to ID. 

We begin with the simplest scenario corresponding to Ji, for which we choose initial 
data Vq such that, in a small neighbourhood of a possible shock at (0, 0), 

(1.16) Vq^ ^ +1 and Vq ^ 0. 
This locally yields the equation 

(1.17) Ut = Ue^oUris) or In (^) = /g f/r(s), where we assume ^ > 0. 
Differentiating this equation in t and, for convenience, renaming the variables 

(1.18) r (— > x and U{r,t) u{x,t) 

yields the necessary equation for such shock wave and blowing up formation: 

(1.19) (t), = ^ or 



UUtt - [UtY = UUxUt. 



Written in the form with the evolution time-variable x and t being the space variable, 

(1-20) ^x = ih*-^]' 

the equation becomes 

(1.21) of parabolic type in the domain of time-monotonicity {ut > 0} 

and parabolic backward in time otherwise. In what follows, we will need carefully control 
the parabolicity condition (11.211) in order to guarantee the well-posedness of our shock 
and blow-up phenomena. In the original independent variables {x,t}, which assume two 
initial conditions at t = 0, fll.lQp exhibits different evolution properties. Note that, in 
terms of t as the time variable, this equation 



;i.22) utt 



is in the normal form in the non-singular set {u 7^ 0}, so it obeys there the Cauchy- 
Kovalevskaya Theorem [21l p. 387]. Hence, for analytic initial data u{x,0) 7^ and 
Ut{x,0), there exists a unique local in time analytic solution u{x,t). On the other hand, 
any existing solution u{x,t), which is analytic in {u 7^ 0}, has a unique local analytic 
continuation at any point, where u ^ 0. 
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Let us make an important comment concerning the type of this new PDE. Written as 
a first-order system for U = {u, v)'^, with v = Uf, f ll.l9p takes the form 



(1.23) Ut = A{U)U^ + F{U), where A{U) 





V 



and F{U) 



It follows that the characteristic equation for the matrix A{U) has the trivial form 

(1.24) det {A{U) - Xh) = =^ = 0. 
Therefore the eigenvalues always coincide, X± =0, so that 

(1.25) (11.231) is not a strictly hyperbolic system, 

and indeed it is uniformly degenerate. In view of (11.251) . application of classic theory 
of hyperbolic systems (see [21 IS]), which recently got a fundamental progress, is not 
possible. In any case, it is worth recalling that formation of singularities for (I1.23P or (11.191) 
and possible extensions of blow-up solutions beyond singularities remain the key (and 
absolutely unavoidable for any difficult nonlinear PDE model) question to be addressed 
in what follows. 

Regardless such nice local regularity and analyticity properties of the PDE (11.191) . our 
main goal is to show that (11.191) admits rather standard "quasi-stationary" self-similar 
shock waves (Section[2]) and other blow-up singularities (Section[3]). In SectionHJ following 
the strategy in [9J, we develop a local theory of formation of shock waves from continuous 
data and show that the PDE (11.191) admits a unique similarity continuation after blow-up. 
In other words, this implies that uniqueness is not lost locally at the singularity blow-up 
points, and this links (I1.19P with the classic model (II. 7p . However, a proper definition 
of entropy solutions to create a global existence-uniqueness theory for (I1.19P remains a 
difficult open problem (it is not clear if such one can be developed in principle). 

1.3. On blow-up similarity solutions with swirl. Let us for a while return to the 
full boxed system in (11.120 . which is a difficult object, so we truly and desperately need 
a simplified model. For instance, in general, (11.120 admits complicated singular solutions 
with blow-up angular swirl. For the NSEs and EEs, such a mechanism was proposed in [7], 
so we omit details, and prescribe those similarity solutions as follows: for fixed a G (0, 1) 
and (T e M, 

Uir,ip,t) = (-ty-f{y,f,), Vir,ip,t) = {-t)-\j{y , ^i) , 

y = '(~tW' = /i = V5 + crln(-t), 

where we introduce blow-up rotation in the angular direction being a logarithmic travelling 
wave. Substituting (I1.26P into (11.121) yields the following first-order system: 



1.27) 



9 + Pg'yV - (^g'f, = gfy 

This is a difficult system, which we do not intent to study it here, though indeed some its 
features will be detected for a reduced simpler PDE model in (I1.19p . 



2. Shock waves by blow-up similarity solutions 

2.1. Existence of self-similar shock formation. Similar to (11. Sp . we look for blow-up 
similarity solutions of (11.191) of the form 

(2.1) u4x,t) = fiy), y = j^y where {1!^)' + t = iff . 
The ODE in (12. ip can be easily reduced to the following one: 

(2.2) yff" = {fY{y + f)-2ff', /(±oo) = ±i, 

which is degenerate at y = and at any zero / = of solutions. Note that, unlike 
(II. 8p . the conditions at infinity have been changed, so that (12. ip is supposed to deliver 
the opposite shock: 

(2.3) \imt-,o- u_{x,t) = S'+(z) = signx. 
In view of the obvious symmetry of (12. 2p : 



(2.4) 



y ^ ~y, 



it suffices to solve the following problem: 

(2.5) yff = (fy(y + f)-2ff for y>0; /(O) = 0, /(+oo) = 1. 
Recall that according to (I1.17P we have to check the sign of 

(2.6) ^ = i-ty^ff > 

which is true if f{y) > and f'{y) > for all > 0. Here, (12. 5p is a pretty standard 
second-order ODE, so we state the final result: 

Theorem 2.1. The problem (12. 5p admits a unique C^([0, oo)) strictly monotone solution. 

Proof. It consists of a few steps, (i) Scaling invariance. One can see that (12. 5p admits a 
group of scalings: if fi{y) is a solution, 

(2.7) then /«(?/) = is a solution for any 0. 

(ii) Finite regularity at y = 0. This can be seen by (still formal) expansion near the origin 

(2.8) f{y) = y + e{y), where e{y) = o{y) as y ^ 0. 

Substituting (12.80 into the ODE fl2.5p and "linearizing" yields the following Euler's ODE: 

(2.9) y^e" -2ye' + e = =^ E{y) = y"'^, where m± = 

Since we are looking for sufficiently smooth and at least classical solutions of the PDE 
(ll.igp . we have to choose the maximal exponent 

(2.10) m = m+ = ^ = 2.61803... G (2, 3) (m_ = 0.38197... < 1). 

The justification of the expansion (12. 8p . (12.90 is performed in a standard manner by 
reducing the ODE to the equivalent integral equations for y = y{f ) (this deletes the trivial 
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Figure 1 . Various solutions of the ODE ()2.5|) with the expansion ()2.1ip for 
different values of yl > 0. 

solution / = 0) and applying Banach's Contraction Principle with a suitable weighted 
functional setting. 

Thus, (I2.5P admits a ID bundle of solutions close to the origin: 

(2.11) f^y) = y + Ay^ + o{y"') for y > 0, 

where A G M is arbitrary. For A = 0, we have the exact unbounded solution 

(2.12) fo{y) = y. 

It follows from the symmetry fl2.4p that the extension of the shock profiles for y < is: 

(2.13) fiy) = y + A\yr~'y + oi\yr) for y<0. 

The parameter A < makes it possible to shoot the necessary constant +1 at infinity, 
which is rather standard and we omit details. □ 

Actually, some necessary properties of solutions of (12.51) can be easily checked numer- 
ically by using a standard MatLab solvers such as the ode45, and this diminishes the 
necessity of a deeper mathematical study. In Figure [T|, we show the monotone depen- 
dence on v4 < and A > of solutions of (12.51) with the expansion (12.111) . For A > 0, 
the local solutions f{y) blow-up in finite yo > that is also easily proved. Note that 
for A < 0, according to (11.211) . the behaviour in {y > 0} of similarity solutions (12.11) are 
locally governed by a well-posed parabolic fiow with a blow-up formation of a shock wave 
in finite time. 

The next Figure [2] demonstrates a general view of similarity shock profiles for both 
y > and y < 0. By the boldface line we denote the unique solution of the problem (12. 5p . 
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Shock wave profiles 
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Figure 2. A general view of similarity shock profiles, 
with the following values in the expansion ( 12.11^ : 

(2.14) A, = -0.96.... 

In addition, in Figure [3] we present the results of further analysis of shock wave profiles 
as (generalized) solutions of the ODE (12. 5p . Namely, we now shoot from y = —50 with 
the following asymptotic expansion of the solutions about the equilibrium / = 1: 

(2.15) f{y) = l + f + ... as y ^ -oo. 

We now use the ode45 solver with the fully enhanced relative and absolute tolerances 

(2.16) Tols = 10-^1 

For i? < 0, we observe that the solutions pass though the singular point y = and 
converge to another larger equilibrium, 

(2.17) f{y)^f+ as y^+oo, where /+ = /+(fi) > 1. 

Note that possible singularities at y = are integrable since even in the worst case the 
solution is bounded, 

=^ f{y) ^ C + y\n\y\ + ... as y 0, 

or the regularity is better, when the right-hand side in (12. 5p vanishes at y = (this case 
actually occurs to give good solutions). 

A different phenomenon is observed for B > 0, where the shock profiles seem have 
finite interfaces at some y = yo < 0. However, next Figure H] explains that this is wrong: 
thought (a) with the enlargement 10~^ continues to convince existence of negative finite 
interfaces, (b) given in the logarithmic scale confirms that the solutions have exponential 
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decay up to ?/ = 0^. Some lines in (b) are disjoint for profiles \f{y)\ ~ 10^^'', which 
suddenly change sign due to numerical errors (recall the guaranteed accuracy fl2.16p ). 
Namely, a more careful look at the equation (12.51) shows that it admits solutions with the 
following typical non-analytic behaviour: 

(2.18) /(y) ~ e"(^ as y ^ 0' , 

where a > is arbitrary due to the scaling symmetry (12.71) . 

Such solutions admit a natural trivial extension by /(y) = for ?/ > 0, thus forming 
C°° shock wave solutions with finite interfaces. Taking such a profile f{y), instead of 
(II. 9p . we get in the limit the reflected Heaviside function: 

(2.19) Um..(x,*) = J/(-x) = {° l^l' 

It is principal to note that, due to the regularity (I2.18P at the interface, the similarity 
solution (12. ip is a C°° function for t < 0. In other words, (I2.19P describes formation of a 
shock wave as t ^ 0~ from classic C°° solutions. Such an evolution process of formation 
of shocks can be taken as a core property for introducing a test on the so-called 6 -entropy 
solutions, i.e., those which are obtained via smooth deformation of shocks appeared. See 
[SI El [ID] , where such an approach is developed for nonlinear dispersion equations such as 

(2.20) ut = {uux)xx and ut = -{uux)xxxx- 

Classic concepts of entropy do not apply to (12.201) . so other approaches are needed. On the 
other hand, a detailed pointwise analysis of formation of gradient blow-up for PDEs such 
as (I2.20p shows that uniqueness is violated there, so a consistent entropy theory is most 
plausibly non-existent in principle; see [9]. Deeper evolution and entropy-like properties 
of such finite interface solutions are out of the scope of this paper. Let us mention that 
the questions on weak compactly supported solutions are in the focus of modern research 
in the EEs area; see [U p. 418]. 

2.2. Existence of rarefaction waves. The corresponding rarefaction waves describe 
collapse of initial singularities. The similarity ones then are defined for alH > 0: 

(2.21) u+{x,t) = F{y), y = f, yFF" = {F')\y-F)-2FF\ F(±oo) = ^l- 

It follows that, in comparison with the shock profiles Fiy) constructed above, the rar- 
efaction ones are given by the refiection: 

(2.22) - F{y) is a solution of the rarefaction ODE ^Tn^ . 

Therefore, the similarity solution in (12.210 is formed by the initial singularity S^{x) = 
— signx and describes its evolution collapse for t > 0. 

Similarly, for data u{x,0) = —H{—x), we obtain similarity solution (12.210 describing 
its collapse. To get a unique solution, two initial functions should be prescribed. It is 
easy to see that ut{x,0) is not a distribution. Indeed, we have, according to (12.210 and 
asymptotics (12.150 . (I2.18p . in the sense of distributions, it is formally valid that 

(2.23) ut{x, t) = -i F'{y)y ao5{x) as t^0+, where Oq = - J F'{y)y dy > 0. 
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Non-symmetric profiles for B>0 and B<0 




y 

Figure 3 . Examples of non-symmetric similarity shock profiles satisfying (I2.15P 
for B < and B > 0. 




However, in view of the asymptotics fl2.15p . the integral in fl2.23p diverges, so that Ut{x, 0) 
is not a bounded positive measure. In the sense of distributions, such a behaviour can be 
characterized by a logarithmic divergence as follows: 

(2.24) utix,t) = F'{y)y r-^ B\\nt\S{x) as t^0+, so ^(x, 0) is a measure. 
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3. Unbounded shocks by similarity blow-up and rarefaction profiles 

More complicated shock-type structure occurs if in (12 .ip we introduce a parameter 
a G (0, 1) looking for similarity solutions (cf. (11.261) ) 

(3.1) u_{x,t) = {-ty''f{y), y = j^ for t < 0, where (3 = 1 -a, 
and f{y) then solves the following more complicated ODE: 

(3.2) /3V/r = /3?/(/0'(/^y + /)+/5(«- 2)2///' + «/'(/'-!) in R. 

Using the same symmetry (12.41) . we again pose the anti-symmetry condition at the origin 

(3.3) /(0)=0 =^ f(-y)^-f(y). 

In view of (I1.17p . we have to require that 

(3.4) f = (_t)-"-i£^£^ > 0, 

which for (13. 2p can be also checked by a Maximum Principle approach. One can see that 
this means that > for a; > implying a monotone in time growth of blow-up solutions 
(a typical feature of parabolic equations, see (ll.2ip ). By a similar shooting technique, we 
prove the following analogy of Theorem 12.11 which includes all necessary asymptotics of 
blow-up profiles: 

Theorem 3.1. The problem (13. 2p . (13.30 admits a family of solutions f{y) > for y > 
satisfying for any A < 0: 



(3.5) fiy)=y + Ay-^ + ... as , 0^ = ^J^L^^^^W^^W^^. 

(3.6) f{y) = Coy'^ + ... ^0+ as y^ +oo, where Cq = Co{A) > 0. 
Let us note that the linearization (12. 8p now leads to the equation 

(3.7) (1 - a)2y V - {a^ - 2a + 2)ye' + (1 - afe = 0, 

whence the same solutions with the exponents m± given in (13. 5p . 

It follows that (13. ip is an unbounded hlow-np solution that forms as t ^ 0~ the following 
singular final time profile (cf. (12. 3p ): 

(3.8) U4..0-) = {f^^^p^^ l>^°^ (._r,0-)=±oo). 

The asymptotics at infinity (13. 6p is also a standard one, which is approved as in reaction- 
diffusion theory; see e.g. [191 Ch. 4] and references therein. 

As usual, we complete our analysis by numerics. In Figure 0, we show a typical blow-up 
profile for a = ^, while Figure E] explains deformation of such profiles f{y) with a. 
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Blow-up similarity profile for a = 0.5 



f(y)=C y-'+.. 




f(y)=y+A y'"++..., A=-1 



Figure 5. A solutions of (fO) . (f3?3]) for a = i and A 
Blow-up similarity profiles for a=0.2, 0.3, 0.4, 0.5 




Figure 6. Solutions of ([321), (03]) for a = 0.2, 0.3, 0.4, and 0.5; A = -1. 



3.1. Rarefaction waves from singular initial data. Similar to (12.211) . the global 
similarity solutions of (11.191) . 



(3.9) 



u+{x,t) =t-"F{y), y = ^ for t > 0, p = 1 - a, 
(j^y'^FF" = /3y{F'f{(3y - F) + (3{a - 2)yFF' + aF\-F' - 1) in 
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corresponding to the reflection F ^ —F in the blow-up ODE (13.21) . describe collapse of 
singular initial data (13. 8p posed at t = 0. We then again arrive at a difficult question on the 
correct entropy-like choice of proper solutions, which reveal some remnants of those ones 
for the EEs (11.11) and hyperbolic systems such as (11.31) . Hopefully, an evolution entropy 
test then can be developed in lines of smooth ^-deformations as in [H [H [10]. However, 
there are examples of higher-order PDEs, for which any entropy-like unique extension of 
a solution after singularity is principally impossible (there are infinitely many extensions 
that are all equally allowed), [9j. We do not know whether such negative trends can be 
oriented to the EEs (fLT]) . 

4. Uniqueness of a local continuation after gradient blow-up 

We now need to consider the principal question on a generic (self-similar) mechanism 
of formation of shocks for the model (ll.lQp . To this end, following [9], we will use the 
same similarity solutions (13. ip . but now with 

(4.1) a < =^ (3 = l-a = l + \a\>l, so that ^-f = < i. 

It can be shown that such solutions (13. ip create continuous "initial data" (cf. (13. 8p ) 

(4.2) M_(x,0") = 



CqX^+I"! for X > 0, 
-Co|a;| 1+1"! for a; < 0, 



where Cq > is some (actually arbitrary by scaling (12. 7p ) constant. Note that (14. 2p is 
still continuous at x = but has a first gradient blow-up there. 

The question is how to extend the data (14. 2 p for t > 0. This is naturally done by 
using the global similarity solution (13. 9p . Unlike the previous cases, we here assume that 
a shock occurs at y = 0, i.e., we put two conditions there 

(4.3) F(0) = Fo > and F'(0) = F^. 

By a rather standard local analysis of the ODE (13. 9p at |/ = 0, it is not difficult to conclude 
that the Cauchy problem (14.30 has a unique smooth solution in the only case 

(4.4) F'{0) = F, = -1 =^ F{y) = Fo-y + S-, + -- 

For other values of F'{0) ^ —1, the solution is not smooth and is singular at y = 0. All 
such singular bundles can be also derived. Therefore, the behaviour of solution F{y) as 
y — > 0^ is very unstable and, even numerically, its construction is not straightforward. 

Before stating the main result, we present in Figure [7] such global similarity profiles 
F(y) for various values of a G [—1,-0.1]. Since the non-monotonicity of F{y) is not 
visible here at all, in the next Figure [8], we show typical shape of such profiles in the 
Cauchy problem posed foTy>S = 0.1, with F{S) = 1 and various F'{6) < 0. Note that 
the absolute minimum points at some y = yo > perfectly corresponds to the Maximum 
Principle for the ODE (13. 9p . Indeed, one has: 

(4.5) F\yo) = 0, Fiyo) > =^ (3V,Fiyo)F" (yo) = \a\F\yo) > 0. 
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y 

Figure 7. Solutions of the ODE dHj}, ([OI) for a = -0.1, -0.2,... -1. 




y 

Figure 8. Typical non-monotone solutions of the ODE dSJD, for a = -0.5. 

By checking the asymptotic properties of solutions of (13. 9p as ?y — +oo, it is not difficult 
to see that all the profiles with the expansion fl4.4l) , which do not vanish by the Maximum 
Principle shown in (14.51) . exhibit the desired behaviour at infinity: 

(4.6) F{y) = Cy^\{l + o{l)) as y ^ +oo. 



Here C is a positive constant but not necessarily the one Cq obtained in (14. 2p by the 
blow-up limit. However, the transition C ^ Cq is uniquely done by the scaling invariance 
( \2.7L so we arrive at: 

Proposition 4.1. For given initial data ( 14.2p . there exists a unique similarity profile F{y) 
such that the global similarity solution fl3.9p satisfies 

(4.7) M_(x,0") = M+(x,0+) m M. 

We call the corresponding profiles {/, F} a similarity global extension pair. Thus, this 
pair is always unique, so the self-similar blow-up solution u_{x,t) for t < has a unique 
extension via the similarity one u+{x,t) for t > 0. 

This justifies a kind of "micro-local" uniqueness theory of shock wave solutions for the 
PDE f ll.lQp . which sounds rather optimistic. Existence, or even a possibility of construc- 
tion, of a global "entropy-like" theory for (11.191) (say, in Oleinik-Kruzhkov-Lax-Glimm- 
Bressan-etc. sense) is quite questionable still. In any case, the above "micro-uniqueness" 
result is rather inspiring bearing in mind that this could affect some similar extension 
properties of the EEs. 



5. Final conclusions: towards the full model of (r, 9?)-shock and blow-up 

FORMATION 

The integral equation (ll.lSp assumes another mechanism of shock formation for data: 
(5.1) Vq^ ^ and Vq ^ 1. 

Then the second term dominates and, on two differentiations, this leads to the PDE: 

W-^) '-'rip - U_/Ut\ _ jj ■ 

Ut\u)t 

This equation is supposed to describe more complicated phenomena of shock and more 
general blow-up formation, where both independent variables (r, Lp) are essentially and 
equally involved. Of course, (15. 2p is uncomparably more difficult than (ll.lQp . and seems 
can give a deeper insight into weak and strong singularity formation for the EEs (and 
also for (11.31) ) in the actual 2D {r, (/j}-geometry. We do not study this equation here. 
Nevertheless, it should be mentioned that equations such as (15. 2p and others must be 
carefully investigated and understood before any serious attack of singularity phenomena 
for Euler equations in 2D and, especially, 3D, which will indeed require other hierarchy 
of some reduced PDE models to introduce and study. 
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